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TRACKING  CLOSELY  SPACED  MULTIPLE  SOURCES 
VIA  SPECTRAL-ESTIMATION  TECHNIQUES 


INTRODUCTION 

The  tracking  of  closely  spaced  multiple  coherent  spatial  sources  is  a  problem  area  in  the  field  of 
radar  where  current  spectral  estimation  techniques  may  find  some  application.  Interest  is  stimulated  by 
two  important  performance  characteristics  reported  for  these  techniques  [1]: 

•  "Superresolution"  which  implies  the  resolution  of  two  or  more  sources  within  the  conven¬ 
tional  beamwidth  of  an  RF  antenna  aperture;  and 

•  "Absence  of  sidelobes"  which  implies  the  resolution  of  two  or  more  sources  of  unequal 
strengths  when  those  sources  are  spaced  more  than  a  beamwidth  apart. 

Both  of  these  advantages  accrue  because  the  Fourier-transform  convolution  of  the  conventional 
antenna-aperture  window  function  is  avoided;  it  is  replaced  by  a  spatial-filter-function  multiplication  [2]. 
The  spatial  filter  is  data  adaptive  to  the  source  scenario,  using  the  available  degrees  of  freedom  of  the 
system  in  an  optimal  LMS  (least-mean-square)  error  sense.  Both  characteristics  are  restricted  to  those 
applications  where  the  processing  system  has  sufficient  degrees  of  freedom  to  model  the  source  scenario 
and  the  sources  have  sufficient  signal-to-noise  ratio  (SNR)  to  be  distinguished.  Also,  both  terms 
should  be  used  with  caution  to  avoid  misinterpretation  with  respect  to  the  conventional  radar  system 
definitions  of  resolution  and  sidelobes,  thus  the  quotation  marks  are  used. 

If  the  RF  sources  in  the  far  field  are  truly  noncoherent,  that  is,  if  their  signals  have  a  zero  correla¬ 
tion  coefficient  at  the  antenna  aperture,  then  impressive  resolution  performance  is  obtainable  via  many 
of  the  optimal  estimation  algorithms  and  techniques  described  in  the  literature.  These  would  include 
the  maximum  entropy  method  (MEM)  of  J.  P.  Burg  (2),  the  maximum  likelihood  method  (MLM)  of  J. 
Capon  and  R.  T.  Lacoss  (21,  the  Wiener  sample  covariance  matrix  inverse  (SMI)  filter  and  its  several 
equivalent  techniques  (1-41,  the  autoregressive  (AR)  least  squares  methods  of  A.  H.  Nuttall  and  L. 
Marple  (5,61,  the  eigenvalue/eigenvector  analysis  techniques  of  R.  L.  Schmidt  and  of  E.  Bienvenu  and 
L.  Kopp  (7,81,  and  others  (91. 

However,  if  the  RF  sources  in  the  far  field  happen  to  be  coherent,  as  is  often  the  case  in  radar 
situations,  then  the  resolution  performance  vs  SNR  for  the  above  algorithms  and  techniques  will  be 
degraded,  often  to  the  point  of  complete  failure  of  resolution.  This  degradation  in  performance  has 
been  addressed  by  several  researchers  (10-12],  and  some  particularly  interesting  and  constructive  results 
have  recently  been  reported  by  Massachusetts  Institute  of  Technology's  Lincoln  Laboratory  [131.  This 
report  reviews  briefly  the  nature  of  the  coherent-source  problem,  discusses  several  of  the  useful  current 
techniques,  and  applies  the  techniques  to  multiple-source  situations  which  include  both  coherence  and 
unequal  strengths  in  order  to  help  illustrate  their  potential  advantages.  All  results  in  this  report  are 
derived  from  computer  simulations  based  upon  the  digital  signal  samples  described  in  the  appendix. 
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The  first  three  terms  are  identical  to  those  in  the  noncoherent  case  and,  taken  alone,  they  are  stationary 
in  space,  because  the  only  space  dependency  involved  is  the  separation  distance  (x,  —  x*).  They  are 
equivalent  to  the  time-domain  situation  of  two  sinusoids  plus  noise,  that  is,  the  spatial  angle  sin(0)  is 
the  dual  of  frequency,  and  the  element  location  is  the  dual  of  time-sampling  instants.  However,  the 
product  term  [the  fourth  term  in  Eq.  (I)]  has  an  amplitude  cosine  variation  that  is  dependent  on  the 
sum  of  the  distances  (x/  -I-  x*)  and  thus  will  not  be  stationary  in  space.  There  is  no  equivalent  in  the 
time  domain.  The  addition  of  the  three  exponent  terms  for  /  -  results  in  the  aperture  amplitude  Junc¬ 
tion  of  Fig.  1.  If  the  coherent  sources  maintain  their  fixed-phase  relationship,  and  if  the  eight-element 
sampling  aperture  does  not  move,  then  Eq.  (1)  will  cause  the  resulting  sample  covariance  matrix  to  be 
non-Toeplitz,  that  is,  the  matrix-element  values  along  diagonals  will  differ  in  accordance  with  the  cosine 
amplitude  function.  Furthermore,  one  will  find  that  this  matrix  has  only  one  unique  eigenvalue,  that 
is,  the  two  sources  are  blended  into  a  single  unique  eigenvector  instead  of  the  usual  two.  Another  way 
of  looking  at  it  is  that  the  two  coherent  sources  are  blended  into  a  single  "virtual  source,"  which  can  be 
quite  different  from  the  point-source/plane-wavefront  model  that  corresponds  to  an  equivalent  time- 
domain  sinusoid. 

These  characteristics  will  usually  result  in  very  unsatisfactory  performance  from  the  processing 
algorithms  listed  above.  A  typical  example  is  shown  in  Fig.  2,  where  we  have  two  30-dB  coherent 
sources  located  at  16°  and  24°  (8°  spatial  separation),  a  fixed  coherent  phase  difference  of  <6c2/i  “  90°, 
a  number  of  snapshots  N  =  1024,  and  an  eight-element  linear-array  sampling  aperture.  The  beamwidth 
of  our  eight-element  array  is  about  16°,  so  that  the  sources  are  separated  by  about  a  half-beamwidth. 
The  algorithm  employed  in  this  example  is  referred  to  as  the  SLC  algorithm,  which  is  described  in  the 
appendix.  Other  examples  may  be  found  in  Ref.  11.  Note  that  the  performance  of  the  SLC  algorithm 
in  Fig.  2b  is  no  better  than  a  conventional  beam  which  is  scanned  through  the  sources.  Fig.  2a.  At  a 
value  of  “  180°  (not  shown),  the  SLC  algorithm  will  indicate  that  two  sources  are  present,  but 
the  locations  are  incorrect. 


(a)  Conventional  beam  scan  (b)  SLC  algorithm 

Fig.  2  —  Spatial  spectrum  estimates  for  30-dB  coherent  sources,  90°  phase,  at  16°  and  24°  (8°  apart), 
8-element  linear  array.  1024  snapshots  covariance  matrix 
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Solutions  to  these  difficulties  may  be  found  in  those  techniques  which  cause  the  product  term  in 
Eq.  (1)  to  average  near  zero,  so  that  the  sample  covariance  matrix  approaches  the  Toeplitz  condition. 
For  example,  consider  a  sideways  movement  of  the  eight-element  sampling  aperture,  such  that  the 
aperture  line  remains  parallel  for  all  snapshot  positions.  Note  that,  whenever  we  shift  the  position  side¬ 
ways  between  snapshots,  we  are  changing  the  received  phase  difference  between  the  two  coherent 
sources,  because  the  distances  from  the  sources  to  the  aperture  phase  center  have  changed.  If  we 
denote  the  incremental  phase  change  per  snapshot  as  A(^„,  then  the  coherent  phase  difference  for  the 
nth  snapshot  referenced  to  the  aperture  phase  center  may  be  written 

=  <l>cO+  ^  (2) 

n— 1 

where  the  starting  phase  difference.  Thus,  as  we  move  sideways  changes  with  time,  and  each 
element  samples  the  aperture-modulation  function  at  different  points  as  it  moves  along.  After  the  sam¬ 
pling  of  one  or  more  complete  cycles  of  the  modulation  function,  similar  interelement-spacing  covari¬ 
ance  terms  become  approximately  equal,  the  sample  covariance  matrix  approaches  the  Toeplitz  condi¬ 
tion,  and  the  resolution  performance  of  the  various  algorithms  becomes  similar  to  their  performance 
against  noncoherent  sources. 

Some  further  thought  on  the  time  rate  of  change  of  in  Eq.  (2)  leads  us  quickly  to  the  conclu¬ 
sion  that  we  have  formed  a  coherent  doppler  shift  between  the  two  sources  via  movement  of  the  phase 
center  of  our  sampling  aperture.  Furthermore,  it  is  obvious  that  the  identical  result  could  be  produced 
if  we  kept  the  position  of  the  sampling  aperture  fixed  but  let  the  corresponding  coherent  doppler  shift 
exist  between  the  two  sources.  Note  that,  for  Fig.  I,  this  would  cause  the  aperture  modulation  function 
to  move  sideways  with  respect  to  the  fixed  sampling  aperture.  Thus,  a  doppler  shift  between  coherent 
sources  and  targets  gives  us  a  second  possibility  for  resolving  them  within  a  beamwidth.  As  indicated 
above,  it  is  desirable  to  have  enough  snapshots  to  sample  one  or  more  complete  doppler  cycles  so  that 
the  covariance  matrix  may  approach  the  Toeplitz  condition.  Figure  3  illustrates  the  easy  resolution  of 
the  same  sources  used  in  Fig.  2,  except  that  source  #2  is  now  given  an  incremental  phase  shift  of  S<t> 
=  1.1°  per  snapshot,  so  that  it  traverses  3.2  doppler  cycles  in  the  1024  snapshots  averaged  for  this  par¬ 
ticular  case.  The  improvement  is  dramatic,  and  performance  is  fully  equal  to  that  in  noncoherent  cases. 
Experience  has  shown  that  performance  is  good  whenever  the  doppler  cyles  exceed  a  value  of  about 
0.75,  with  adequate  sampling  of  the  aperture-modulation  envelope  function. 

FORWARD-BACKWARD  SCBAPERTURE-SHIFT  SOLUTION 

If  a  doppler-cycle  shift  is  not  available,  then  we  are  forced  to  address  the  difficult  fixed-phase 
coherent  case.  The  best  technique  found  to  date  is  to  form  a  reduced-dimension  sample  covariance 
matrix  via  synthetic  shifting  of  a  reduced-dimension  sampling  subaperture  along  the  elements  of  the 
array.  If  L  represents  the  number  of  elements  in  the  sampling  subaperture,  and  K  is  the  total  number 
of  elements  in  our  linear  array  antenna,  then  we  can  form  a  reduced  L  x  L  sample  covariance  matrix 
averaged  over  (K  -  L  +  1)  samplings  per  array  snapshot.  Such  synthetic  movement  of  a  subaperture 
is  essentially  a  fractional  doppler-cycle  shifting  of  the  source  phasing  at  the  phase  center  of  the  subaper¬ 
ture  as  it  moves  along  the  array.  Furthermore,  the  averaging  is  usually  improved  significantly  if  our 
sampling  subaperture  is  reversed  and  moves  across  the  array  in  the  opposite  direction.  This  double 
synthetic  sampling  movement  is  illustrated  in  Fig.  4,  where  the  reader  will  recognize  its  similarity  to  the 
forward-backward  linear  prediction  techniques  described  in  the  literature  15,6).  Note  that  the  resulting 
reduced  L  x  L  sample  covariance  matrix  will  be  symmetric  but  still  non-Toeplitz.  Figure  5  illustrates 
the  resolution  of  two  30-dB  coherent  sources  spaced  4°  apart  for  L  =  8  and  A  —  16  elements,  utilizing 
the  SLC  algorithm  for  various  coherent  fixed  phases  of  0,  10,  20,  90,  160,  170,  and  180°. 
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Fig.  3  —  Same  coherent-source  situation  as  in  Fig. 
2,  except  that  doppler  shift  of  1.1°  per  snapshot  is 
inserted 
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FORWARD  SHIFT  BACKWARD  SHIFT 

Fig.  4  —  Illustration  of  forward-backward  shift  movement 
for  a  reduced-dimension  sampling  subaperture  where  Z.  =  4, 
along  a  linear  antenna  array  o(  K  15  elements 


Although  it  is  not  apparent  in  Fig,  5,  there  is  a  substantial  performance  dependency  upon  the 
fixed-phase  difference  between  the  two  coherent  sources.  When  (t>(2n  ""  ±90®,  performance  is  about 
equal  to  that  in  the  noncoherent  case,  because  the  forward-backward  averaging  effect  is  maximum. 
However,  when  <t>c2n  —  0°  ot  180°,  then  the  forward-backward  averaging  effect  is  minimum,  because  of 
the  symmetry  of  the  aperture-modulation  envelope  function,  and  performance  suffers  with  respect  to 
that  in  the  noncoherent  case. 

The  forward-backward  subaperture-shift  solution  can  also  be  applied  to  the  single-snapshot  case. 
Note  that,  when  only  one  snapshot  of  data  is  available,  we  have  a  fixed-phase  coherent  case  even  if  the 
sources  are  nominally  noncoherent.  The  single-snapshot  case  is  the  most  difficult  of  all,  because  the 
particular  noise  voltages  frozen  into  that  snapshot  may  have  a  large  variance,  bias  offsets,  or  spurious 
peaks.  To  reduce  the  effects  of  the  noise,  it  is  always  desirable  to  have  a  number  of  snapshots  available 
for  processing. 

EIGENVALUE/EIGENVECTOR  ALGORITHMS 

When  a  signal  is  known  to  consist  of  pure  sinusoids  in  white  noise,  an  appropriate  procedure  to 
find  the  unknown  frequencies  and  powers  of  the  sinusoids  in  the  signal  is  the  Pisarenko  spectral- 
decomposition  procedure  [14].  Much  of  the  current  interest  in  eigenvector  methods  can  be  traced  to 
Pisarenko's  work.  Recent  formulations  specifically  tailored  to  the  spatial  direction-finding  problem  have 
been  described  by  Schmidt  [7]  and  Bienvenu  and  Kopp  [8].  These  techniques  are  based  upon  the 
geometrical  vector  property  that  the  "noise"  eigenvectors  of  a  covariance  matrix  are  orthogonal  to 
source  direction  vectors  and,  therefore,  will  produce  zeros  (nulls)  at  source  locations. 

The  particular  version  used  for  processing  the  simulated  snapshot  data  discussed  in  this  report  can 
be  derived  from  the  constrained  optimized  beam  scan  algorithm  [1]  or  its  identical  counterpart,  the 
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Fig.  5  —  .Spatial  .spectrum  estimates  for  two  30-dB  coherent  sources  located  at  18'  and  22“  (0.5 
beamwidths  apart),  lb-element  linear  airav.  32  snapshots  at  each  ft.xed-phase  value  of  0.  10.  20.  90,  160. 

170,  and  180°  utilizing  subaperture  /.  -■  •' 

MLM  algorithm  of  Capon  and  Lacoss  (21.  where  the  spatial  spectrum  would  be  estimated  by  the  output 
residual  power  P„  from  the  optimized  adapted  array  weights, 

P„  =  W;'MW„.  (3) 

where 

W„  =  optimized  weights  =  S*. 


M  =  covariance  matrix  estimate, 

S*  =  mainbeam  direction  steering  vector,  and 
fi  =  a  scalar  quantity. 

Under  the  zeroth-order  gain  constraint,  we  require  that  S'  =  1,  whereupon 

fi  =  (S'  M-'S*)-'. 

Substituting  /x  and  W„  into  Eq.  (3)  then  results  in 

P  =  - ^ - 

"  S' M-'S*' 


(4) 

(5) 


If  the  steering  vector  S*  for  a  given  covariance-matrix  inverse  is  swept,  P„  will  estimate  the  spatial 
spectrum. 
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From  the  properties  of  the  associated  eigenvector  orthonormal  modal-matrix  transformation, 
defined  as  the  matrix  Q,  we  can  rewrite  the  covariance-matrix  inverse  in  the  form 


M  '  =  Q'[q‘  M  Q'I  O’ 


(6) 


=  Q' 


where  /3,^  is  the  /th  eigenvalue  of  the  covariance  matrix.  This  then  permits  us  to  form  a  simple  sum¬ 
mation  of  eigenvector  power  beams  referenced  to  the  receiver  noise-power  level  /Jq, 


where 


00 

K 

00^ 

Po 

=  z 

k-\ 

0i 

(e). 


(7) 


gi,(9)  =  S'e*. 


(8) 


The  A:th  eigenvector  beam,  g^iO),  is  defined  by  the  <rth  eigenvector  e*.  If  we  use  all  K  beams,  we  have 
an  identity  with  the  constrained-optimized  algorithm;  but  if  we  select  only  the  noise  eigenvector  beams 
based  upon  a  threshold  criterion  related  to  (3^,  then  we  should  get  excellent  nulls  at  all  source  positions, 
because  the  noise  eigenvectors  are  orthogonal  to  the  source  direction  vectors  (7). 

Under  ideal  simulation  conditions,  such  as  are  assumed  in  this  report,  the  threshold  criterion  is  a 
simple  set  value,  because  the  noise  eigenvalues  are  well  behaved  and  have  a  small  spread  in  magnitude. 
However,  in  practical  situations  where  system  errors  and  imperfections  tend  to  widen  the  spread  in 
noise  eigenvalues,  one  must  develop  a  more  sophisticated  criterion  tailored  to  the  specific  application. 

An  additional  dividend  from  this  algorithm  is  that,  once  the  locations  of  the  sources  are  deter¬ 
mined,  then  their  relative  strengths  may  be  evaluated  from  Eqs.  (5)  or  (7),  that  is,  the  MLM  algo¬ 
rithm.  This  pseudo-linear  superposition  property  is  dependent  upon  the  reciprocal  of  a  matrix  inverse 
being  equal  to  the  sum  of  the  reciprocals  of  the  constituent  matrix  inverses.  Other  power  estimation 
techniques  with  better  characteristics  are  available  [7,13,14].  It  should  be  pointed  out  that  estimation 
of  source  power  is  a  key  issue  in  the  use  of  high-resolution  techniques  to  sort  out  and  reject  false  peaks 
in  the  spectral  estimates.  This  becomes  especially  important  when  only  a  few  snapshots  are  available 
and  the  SNR  is  low. 

RESOLUTION  PERFORMANCE:  MIXED  MULTIPLE  SOURCES 

Multiple  source  situations  which  include  both  coherence  and  unequal  strengths  constitute  an 
excellent  vehicle  for  testing  estimation  algorithms,  and  Fig.  6  illustrates  an  example  in  which  we 
attempt  to  estimate  six  sources  using  an  eight-element  linear-array  sampling  aperture.  This  particular 
test  case  includes  two  10-dB  sources  95%  correlated,  at  -54°  and  -42°  (0.57  beamwidths  apart),  and 
four  noncoherent  sources  of  -9-,  30-,  10-,  and  20-dB  strengths,  located  at  0°,  15°,  22°,  and  30°.  Note 
in  Fig.  6a  that  the  conventional  beam  scan  (Fourier  transform)  is  hopeless,  as  expected,  because  of  its 
beamwidth  and  sidelobe  limitations.  In  addition,  the  simple  SLC  algorithm  has  failed  to  resolve  the 
highly  correlated  sources  and  the  two  weaker  noncoherent  sources  located  at  0°  and  22°.  However,  in 
Fig.  6b  the  eigenanalysis  algorithm  has  located  all  six  sources  correctly  via  the  "noise  eigenvectors"  and, 
in  addition,  permits  us  to  evaluate  the  relative  strengths  of  the  noncoherent  sources.  This  algorithm 
has  consistently  demonstrated  more  than  an  order-of-magnitude  improvement  in  SNR/resolution  per¬ 
formance  over  the  older  SLC  algorithm. 

Figure  6c  illustrates  a  typical  disaster  that  occurs  whenever  one  attempts  to  work  with  a  Toeplitz 
matrix  constraint  (autocorrelation  matrix)  in  a  coherent  source  situation.  Interestingly,  not  only  is  the 
coherent  source  information  distorted,  but  the  noncoherent  source  information  is  badly  affected  as  well. 
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(c)  Eigenvector  algorithm,  Toeplitz  matrix  constraint 

Fig.  6  (Continued)  —  Spatial  spectrum  estimate  for  six  sources:  two  10-dB  sources  95%  correlated  at  -54° 
and  -42°  (0.57  beamwidths  apart);  four  noncoherent  sources  of  -9-,  30-,  10-,  and  20-dB  strengths  located 
at  0°,  15°,  22°,  and  30°;  8-element  linear  array,  sample  covariance  matrix  of  1024  snapshots 

A  MOVING-SOURCE  TRACKING  SIMULATION 

The  second  example  involves  four  moving  sources,  of  unequal  strength,  traveling  in  criss-crossing 
patterns.  The  true  paths  for  each  source  vs  time  are  indicated  in  Fig.  7.  Sources  51,  52,  and  53  are 
coherent  sources  of  low  {— 10-dB)  strength;  51  is  the  reference  source;  52  is  95%  correlated,  with  a 
slow  doppler  shift  of  0.027°  per  snapshot;  and  53  is  95%  correlated,  with  a  moderate  doppler  shift  of 
0.586°  per  snapshot.  Source  Jq  is  a  noncoherent  20-dB  source,  intended  to  represent  a  strong  jammer. 
The  time  scale  of  90  units  represents  59850  snapshots  taken  at  equal  sampling  intervals,  so  that  665 
snapshots  are  sampled  within  each  time  unit  along  the  scale.  Thus,  for  each  time  unit,  source  52 
advances  18°  in  phase  and  source  53  advances  390°  in  phase  with  respect  to  51.  The  sampling  aperture 
used  in  this  example  consists  of  a  12-element  linear  array  with  half-wavelength  spacing. 

A  conventional  Fourier-transform  performance  basis  for  this  example  is  given  in  Fig.  8,  where  we 
have  three  different  output  display  plots  based  upon  scanning  beam  processing.  Figure  8a  illustrates  a 
three-dimensional  (3D)  display  output  obtained  from  a  conventional  scanning  beam  of  uniform  illumi¬ 
nation,  identical  to  a  Bartlett-window  Fourier  transform,  wherein  each  time-unit  plot  cut  is  computed 
from  the  sample  covariance  matrix  averaged  over  its  665  snapshots, 

P„-S*'MS,  (9) 

where  M  is  the  sample  covariance  matrix  and  S  is  the  beamsweep  steering  vector  used  to  generate  the 
display  plot.  As  expected,  this  simple  Fourier  output  is  dominated  by  the  strong  source  Jq  to  such  an 
extent  that  it  is  impossible  to  see  the  weak  sources,  let  alone  track  them.  If  one  attempts  to  utilize  a 
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Fig.  7  —  True  paths  used  in  tracking 
simulations  for  four  moving  sources; 
SI,  SI.  and  S3  are  coherent  and  of 
-10-dB  strength;  is  a  strong  20-dB 
noncoherent  source 


low-sidelobe  window,  as  illustrated  in  Fig.  8b,  the  results  are  no  better,  because  of  the  increase  in 
beamwidth.  Whenever  one  has  a  dominating  strong-source  situation  like  this,  a  good  conventional 
technique  to  employ  is  pattern  subtraction  of  the  strong  source,  and  Fig.  8c  illustrates  the  considerable 
improvement  available.  Note  that  the  paths  of  the  three  weak  sources  are  now  crudely  outlined.  Pat¬ 
tern  subtraction  is  similar  to  the  use  of  a  spatial  filter  which  nulls  out  the  strong  source.  We  need  to 
know  the  direction  angle,  9o,  of  the  strong  source,  and  this  can  be  obtained  via  automatic  tracking  of 
source  Jq  in  our  example.  Knowing  0q,  we  can  form  a  unit  filter  vector  F„,  where  the  fcth  element 
component  may  be  written 

(10) 

where 


«o 


sin  9o. 


The  correlation  coefficient  o„  between  F„  and  the  snapshot  signal-sample  vector  E„  defined  in  the 
appendix  is  given  by 

(11) 

The  filtered  signal  sample  vector  E,J  is  then  obtained  by  subtracting  the  filter  direction  components  in 
Gram-Schmidt  fashion,  that  is, 

E;  =  E„-a,F„.  (12) 

From  E,;  one  can  form  the  filtered  sample  covariance  matrix  M'  and  substitute  into  Eq.  (9)  to  compute 
the  time-unit  plot  cuts  shown  in  Fig.  8c. 
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For  comparison  against  the  Fourier-transform  processing  results  shown  in  Fig.  8,  the  results  com¬ 
puted  from  eigenanalysis  algorithm  processing  of  the  same  snapshot  data  are  illustrated  in  Fig.  9.  In 
Fig.  9a,  we  have  the  3D  display  computed  from  the  full  12-element  sample  covariance  matrix  with  no 
subaperture  shift.  Note  that  the  correct  tracks  of  all  three  weak  sources  are  clearly  indicated,  even  in 
the  presence  of  the  dominating  strong  source  Jq.  In  addition,  superresolution  of  these  weak  sources  is 
evidenced  in  about  30%  of  the  time  display  and  contributes  considerably  to  delineating  the  paths  near 
crossovers. 


Fig.  9  —  Eigenanalysis  algorithm  estimates  for  the  four  moving  sources  referenced  in  Fig.  8;  same  12-element  linear-array 
snapshot  data;  each  lime-unit  plot  cut  is  processed  from  665  snapshots. 
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It  is  of  interest  to  note  that  sources  SI  and  S2  are  not  resolved  for  times  T  =  0,  1,  and  2.  The 
reason  for  this  may  be  seen  if  we  examine  the  coherent  source  information  listed  in  Table  1,  which 
shows  that  these  two  sources  start  out  in  phase  at  a  separation  of  0.5  beamwidths.  There  is  a  noticeable 
peak  modulation  along  the  paths  of  SI  and  S2,  corresponding  to  the  slowly  changing  phase  difference 
between  them.  This  modulation  persists  even  when  they  are  separated  by  several  beamwidths. 


Table  I  —  Coherent  Sources  SI  and  S2 


Time 

T 

Spacing 

(beamwidths) 

Average  Phase 
Difference 
\<f>  (deg) 

0 

0.50 

0 

0.55 

18 

0.60 

36 

0.65 

54 

0.70 

72 

0.75 

90 

0.80 

108 

0.85 

126 

8 

0.90 

144 

mm 

0.95 

162 

■■ 

1.00 

180 

The  3D  display  output  can  usually  be  improved  significantly  for  coherent  sources  by  the  incor¬ 
poration  of  some  forward-backward  subaperture-shift  processing,  as  described  earlier.  Figure  9b  illus¬ 
trates  this  improvement  for  the  current  example,  with  a  subaperture  of  eight  elements  used  to  process 
the  same  12-element  data  snapshots.  Note  that  not  only  are  the  peaks  higher  and  cleaner,  but  the 
modulation  of  the  peak  amplitudes  referred  to  above  has  been  largely  eliminated. 

Thus  far  in  this  final  example,  each  plot  cut  was  processed  from  the  full  665  snapshots  per  time 
unit.  If  we  reduce  the  number  of  snapshots  processed,  then  the  averaging  in  the  sample  covariance 
matrix  is  proportionally  reduced,  the  variance  of  our  estimates  increases,  and  false  alarms  or  spurious 
estimates  will  increase.  In  brief,  the  3D  display  output  will  tend  to  become  noisy.  Figure  10  illustrates 
this  noisy  effect  caused  by  a  reduction  to  only  32  snapshots  processed  per  time  unit.  Note  that  the 
track  of  the  strong  source  Jq  is  not  bothered  much  by  the  reduction,  because  of  its  high  SNR. 

More  dramatic  differences  between  conventional  Fourier  transform  processing  and  eigenanalysis 
processing  could  readily  be  demonstrated  by  an  increase  in  the  strengths  of  sources  S2  and  S3  in  our 
example  or  by  the  the  addition  of  more  sources  of  unequal  strengths. 

CONCLUSIONS 

Current  spectral-estimation  techniques  have  been  applied  to  multiple-source  situations,  which 
include  both  coherence  and  unequal  strengths,  in  cases  designed  to  illustrate  two  of  the  important  per¬ 
formance  characteristics  of  interest  in  radar  or  passive  tracking:  superresolution  and  apparent  absence  of 
sidelobes.  The  simulations  demonstrate  a  resolution  capability  that  is  very  impressive  compared  to  con¬ 
ventional  techniques,  and  they  have  optimistic  implications  regarding  the  simultaneous  accurate  track¬ 
ing  of  closely  spaced,  mixed,  multiple  sources  or  targets.  In  addition  to  these  results,  recent  work 
reported  by  other  researchers  [13]  also  indicates  that  the  level  of  performance  achieved  with  modern 
spectral-estimation  techniques  has  reached  the  point  where  these  techniques  invite  in-depth  considera¬ 
tion  for  real-world  tracking-system  applications. 
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Fig.  10  —  Elgenanalysis  algorithm  estimates 
with  the  same  procc.soing  as  Fig.  9b.  except  that 
the  number  of  snapshots  processed  per  time- 
unit  plot  cut  is  reduced  to  32 


It  is  recognized  that  a  number  of  serious  issues  and  problem  areas  remain  to  be  addressed.  For 
example; 

•  The  basic  reliability  of  source  or  target  location  in  a  statistical  sense  under  practical  track¬ 
ing  conditions.  This  includes  threshold  criteria  and  false  alarm  problems. 

•  The  compatibility  of  mixing  with  conventional  range  and  doppler  signal  processing  to 
achieve  improved  four-dimensional  radar  space-resolution  performance. 

•  The  robustness  of  the  model-dependent  data-adaptive  processing  when  the  source 
scenario  exceeds  the  available  degrees  of  freedom,  causing  saturation  degradation. 

•  The  basic  array  and  channel  error  and  accuracy  considerations.  This  includes  the  often- 
overlooked  requirement  for  calibration  of  the  element  or  beam  patterns  for  all  scan  angles 
and  bandwidths. 

•  The  extensive  impact  upon  overall  system  design,  from  the  antenna  to  the  interactive  out¬ 
put  display.  An  all-digital  system  is  clearly  indicated. 

•  The  cost  effectiveness  compared  to  conventional  techniques. 

Obviously,  considerable  research  effort  is  required  to  achieve  the  desired  goals,  but  the  outlook  is 
optimistic  because  of  the  rate  of  progress  which  has  characterized  these  techniques  in  recent  years. 
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Appendix 

SIGNAL  SAMPLES  AND  SLC  ALGORITHM 


Consider  a  simple  linear  array  of  K  equally  spaced  elements.  The  received-signal  samples  will  be 
correlated  in  both  space  and  time,  giving  rise  to  a  two-dimensional  data  problem,  but  we  convert  this  to 
the  spatial  domain  only  by  assuming  that  narrowband  filtering  precedes  our  spatial-domain  processing. 
Signal  samples  are  defined  as  aperture  snapshots,  that  is,  one  simultaneous  sampling  of  the  aperture  sig¬ 
nals  at  all  array  elements.  The  nth  snapshot  signal  sample  at  the  A:th  element  will  consist  of  indepen¬ 
dent  Gaussian  receiver  noise,  plus  the  two  coherent  source  voltages, 

£■*,  -  -I-  (Al) 

where 

'<'/*»  -  +  <t>J. 

lit  .  « 

M,  “  sm  tf,, 

and  X  =  wavelength,  x*  =  position  of  A:th  element  from  center  in  wavelengths,  =  spatial  location  of 
/th  source,  Aj  -  amplitude  of  /th  source,  =  random-sampled  phase  of  both  sources  for  the  nth 
snapshot,  <t)c\„  =  0  so  that  source  #1  becomes  our  coherent-source  reference,  and  then  incor¬ 
porates  the  coherent  phase  difference  between  the  two  sources.  A  nontrivial  point  is  that  the  coherent 
phase  difference  is  ref^erenced  to  the  center  of  the  array. 

A  convenient  vehicle  for  application  of  the  matrix-inverse  type  of  algorithm  is  the  generic  SLC 
(sidelobe-canceler)  antenna  configuration  illustrated  in  Fig.  Al.  The  reader  will  note  its  resemblance  to 
a  linear  prediction  filter  [1].  Throughout  the  report,  this  configuration  is  referred  to  simply  as  the  SLC 
algorithm.  The  spatial  spectrum  estimates  are  the  reciprocals  of  the  adaptive-filter  patterns,  which,  in 
turn,  are  computed  from  the  optimum  weights,  W^.  We  simply  compute  W„  from  the  inverse  of  the 
sample  covariance  matrix  [3], 


W,  - 

(A2) 

S*'-  (0,  0,  0,  0,  0,  0,  0,  1], 

(A3) 

(A4) 

M.  -  Ie;  e;i, 

(A5) 

where  E„  is  the  /7th  snapshot  signal-sample  vector,  whose  element  components  are  given  by  Eq.  (Al), 
M„  is  the  /7th  snapshot  contribution  to  the  covariance  matrix,  M  is  the  sample  covariance  matrix  aver¬ 
aged  over  snapshots,  S*  is  the  quiescent-weight  steering  vector,  and  ;u  is  a  scalar  quantity.  Note  that, 
in  true  sidelobe-canceler  fashion,  the  steering  vector  S*  injects  a  zero  weight  on  every  element  except 
for  the  end  element,  thus  causing  the  quiescent  pattern  of  the  array  to  be  that  of  the  single  end  ele¬ 
ment.  Several  examples  of  the  use  of  this  SLC  algorithm  against  noncoherent  sources  can  be  found  in 
Ref.  [1). 
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